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We investigate the difference between classical and quantum dynamics of coupled magnètic 
dipoles. We prové that in general the dynamics of the classical interaction Hamiltonian differs 
from the corresponding quantum model, regardless of the initial state. The difference appears as 
non positive-definite diffusion terms in the quantum evolution equation of an appropriate positive 
phase-space probability density. Thus, it is not possible to express the dynamics in terms of a 
convolution of a positive transition probability function and the initial condition as can be done in the 
classical case. We conclude that the dynamics is a quantum element of NMR quantum information 
processing. There are two limits where our quantum evolution coincide with the classical one: the 
short time limit before spin-spin interaction sets in and the long time limit when phase diffusion is 
incorporated. 
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I. INTRODUCTION 

Recent work in quantum information theory has sug- 
gested that quantum computers are more powerful than 
their classical counterparts In quantum communi- 

cation, there already exist algorithms which have been 
proved to out-perform their classical counterparts 
The situation in computation is not as clear: we know 
of problems which have quantum algorithms that are ex- 
ponentially faster but only when compared to the known 
classical ones [|ï],|]|||6||| (not the optimum ones). The 
special power of quantum computers is only a conjecture 
as we have no proofs that we cannot simulate efficiently 
quantum systems using classical computers. Although 
most physicists would believe this efficient simulation to 
be impossible, this is at the foundation of the distinction 
between classical and quantum computation. 

If quantum computers are indeed more powerful than 
classical ones, could we pinpoint the origin of this power 
to one or a few elements in the quantum mechanical the- 
ory? In the 'folklore', the power of quantum compu- 
tation has been attributed to entanglement. In quan- 
tum computers where the initial states are pure it has 
been claimed Jï(],[ïï| that the presence of entanglement 
distinguishes quantum and classical algorithms. Indeed 
the evolution of a quantum system starting in a pure 
state and evolving unitarily without entanglement (which 
could oceur if there were no interactions between the 
components) can be efficiently simulated by a classi- 
cal system; on the other hand, a classical simulation of 
genèric quantum evolution for a pure state has no known 
efficient algorithm. At the basis of this argument is the 
ability to efficiently simulate a system by a classical Com- 
puter. 

The argument employed in Ref. |ïï|] , which usues the 
existence of an efficient classical simulation in the ab- 
sence of entanglement, does not carry through when the 



initial state is mixed, (that is not a pure state). That is, 
for some highly mixed state |ï^] such as the state of nu- 
clear spins present in liquid-state NMR, we do not know 
how to efficiently simulate the evolution of the system on 
a classical computer. For these states, the density ma- 
trix can be represented as a sum of separable states with 
positive coefficient (to be interpreted classically as prob- 
abilities to be in the respective states). However, under 
genèric unitary evolutions, the choice of separable states 
must change. No efficient algorithms exist to relate the 
initial separable states to the final ones for an increasing 
number of spins. 

We do not yet have a genèric quantitative measure for 
entanglement, although we do have a measure for absence 
of entanglement. Pure states are defined as being separa- 
ble, or non-entangled, if they can be expressed as prod- 
uets of subsystem (such as qubit) wave functions. For 
mixed states, this notion is generalized to the existence 
of at least one expansion of the state in terms of separa- 
ble pure states with positive coefficients. Thus an equal 
mixture of the maximally entangled state of two spins 
(Bell states) does not contain any entanglement because 
this density matrix can be re-expressed in terms of sepa- 
rable states (for example the computational states with 
equal probability) . The separable states of spin-half sys- 
tems could be described at a given time as a probability 
distribution of a set of classical tops. 

The notion of entanglement for mixed states has been 
developed in the context of quantum communication. 
One definition corresponds to the number of maximally 
entangled states that can be extracted ensemble of these 
states flÏ3| . But a computation is inherently a dynamical 
process, and we do not know in general how to describe 
the evolution of one mixed separable state to another 
using an efficient classical description, in contradistinc- 
tion to the pure state case. The quantum device can thus 
provide some information more efficiently that a classical 
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device could. 

If mixed states are used as initial states of a quan- 
tum computer, entanglement does not seem to play the 
essential role in distinguishing quantum and classical al- 
gorithms as it might if we used pure states. A particular 
example of a quantum algorithm without a known cffi- 
cient classical counterpart is the one given in Ref. JÏ4) . 
The algorithm gives the distribution of eigenvalues of a 
quantum Hamiltonian. It uses as an input state an ex- 
tremely mixed state, one with a single qubit in a pure 
state (a pseudo-pure state Q would also do), and all 
other qubits are maximally mixed. This algorithm uses 
extremely mixed states but can still outperform known 
classical algorithms. For a small number of qubits there 
will definitely be no entanglement if a pseudo-pure state 
of the first bit is used; as we increase the number of 
bits there may or may not be entanglement. But even 
if it happens at the n th qubit, the algorithm will not go 
through any phase transition; thus, it would be meaning- 
less to refer to the algorithm as classical before the pres- 
ence of entanglement and as quantum afterwards. What 
distinguishes this algorithm from the classical analog is 
that the rules for transforming the density matrix are the 
quantum rules, and we do not know how to efficiently 
simulate them by the classical rules. 

The algorithm in Ref. 0] is especially relevant in the 
context of recent discussions of experiments in quantum 
information processing using liquid-state NMR technol- 
ogy JÏ0yÏ6|] , The algorithm in Ref. [Q could be imple- 
mented in liquid-state NMR. The authors of |ï(| com- 
mented: "The results in this Letter suggest that current 
NMR experiments are not true quantum computations, 
since no entanglement appears in the physical states at 
any stage." This statement makes the assumption that 
entanglement is the necessary element of quantum com- 
putation following the suggestion in ]ï"ï|] . I n the same 
paper [jïo| however it is recognised that it may not be so 
easy to separate quantum dynamics and entanglement 
when trying to pinpoint the power of quantum copu- 
tation; "The results in this Letter suggest that current 
NMR experiments are not true quantum computations, 
since no entanglement appears in the physical states at 
any stage. We stress, however, that we have not proved 
this suggestion, since we would need to analyze the power 
of general unitary operations in their action on separable 
states. To reach a firm conclusion, much more needs to 
be understood about what it means for a computation to 
be a "quantum" computation." However the claim that 
the evolution of unentangled pure states can be efficiently 
simulated by a classical computer JÏÏJ] does not carry 
through to mixed states. The power of quantum com- 
putation can come from properties of the dynamics, not 
the state. This was also recognised by Schack and Caves 
|Ï6f . Indeed, the real origin of the criteria in (ï^] is the 
dynamical evolution of the system not the state itself. If 
we apply the type of unitary transformation used in ]TÏ| ] 
to a highly mixed state (so that entanglement might not 
appear), it is as hard to simulate on a classical computer 



as when the state is initially pure. This point was not 
considered in Ref. |ïü| . 

Recognizing this fact, Schack and Caves attempted to 
explain some liquid-state NMR experiments using classi- 
cal dynamics, without success JT^] . They did not derive 
an equation of motion for the behavior of the spins but 
rather provided a model which described the effect of 
"gates" on the states. Their model did predict an expo- 
nential decay of the signal as a function of the number of 
gates going as (1 + for n the number of qubits 

in the experiment and g the number of gates. Even in 
the three qubit experiments they commented upon, the 
gates defined by Schack and Caves are not the physical 
gates implemented in the experiment as this would have 
ruled out their model. In the seven-qubit experiment of 
Ref. ]ÏSt| , the model predicted a decrease in the signal- 
to-noise ratio of ~ 10 40 smaller than observed: a number 
which unquestionably rules out the model. Discrepàncies 
of their theory with older NMR experiments, although 
not called quantum computing at the time, are plenti- 
ful in the NMR literature, and the readers can review 
experiments performed twenty years ago pp|pl[ |. 

What does make the quantum dynamics so hard to 
simulate? Could there be other classical models which 
explain NMR experiments? Can we understand the ori- 
gin of the discrepancy between classical and quantum 
evolution? In this paper we compare the evolution of 
two coupled spin-half partides under quantum and clas- 
sical evolution. The work of the last fifty years in NMR 
shows the consistency of experimental results with quan- 
tum mechanics and the failure to find a classical descrip- 
tion: "The dynamics of isolated spins can be understood 
in terms of the motion of classical magnetization vectors. 
To describe coupled spins, however, it is necessary to have 
recourse to a quantum mechanical formalism where the 
state of the system is expressed by a state function or, 
more generally, by a density operator." [^2| Here, we will 
give the explícit origin for this difference for the simplest 
choice of a classical model: the one with the same Hamil- 
tonian as the quantum model. In the next section we 
derive the evolution equation for classical and quantum 
interacting spin-half partides, explicitly demonstrating 
that the classical theory (with the same Hamiltonian) 
cannot reproduce the quantum equations. We then dis- 
cuss implications of these equations and draw conclu- 
sions. 



II. CLASSICAL AND QUANTUM DYNAMICS OF 
COUPLED SPIN-HALF PARTICLES 

We investigate the equations of motion of classical and 
quantum spin-half particle and show that the quantum 
behavior is fundamentally different from the classical one. 
We show that even in the cases where the density is highly 
mixed (non-entangled) the evolution leads to different 
observables quantities. There is some ambiguity in ex- 
actly what is meant by the classical dynamics of such a 
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system. We must agree on some ground rules to make 
a meaningful comparison. Semiclassical dynamics has a 
long history and some rules have been established p3[ . 
We will assign a classical analogue for a quantum prob- 
lem by demanding that the same functional form of the 
Hamiltonian be used but with the corresponding clas- 
sical phase space variables substituted for the quantum 
canonical operators. 

In the case of a spin-half system this would appear 
to present some problems, but the situation is clearer 
if we always work in the irreducible representations of 
the total angular momentum of the system. There are 
two subspaces corresponding to total angular momentum 
quantum number s = 1 and s = 0. The dynamics of a 
single spin-half system will conserve angular momentum. 
However, for two spin-half systems, with arbitrary one- 
and two- qubit interactions, angular momentum need not 
be conserved and these subspaces will become coherently 
mixed. It will suffice however to consider one particular 
two qubit gate that does conserve angular momentum, 
which is an important two qubit gate for quantum com- 
puting and for which the classical and quantum dynamics 
are completely different except on a short time scale. Fur- 
thermore this is precisely the two qubit gate accessible in 
NMR quantum computing. 

In order to produce entangled states of single spin-half 
systems, a variety of possible interactions could be used; 
however, in NMR the natural intcraction is of the form 
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with H = 1 and the subscript 2 indicating that the inter- 
action is between two qubits. The dynamics of two spin 
states that follows from this Hamiltonian is given by a 
unitary operator ^(í) where 



U 2 (t) = exp 



4 z z 



(2) 



In addition to this two-spin unitary operator or 'gate', 
genèric quantum computation needs single spin dynam- 
ics. This is easily generated by the scalar coupling of a 
spin-half magnètic dipole with an applied magnètic field. 
The Hamiltonian for these single spin rotations for the 
spin labelled i is 



Hx = B(t) ■ <J 



(i) 



(3) 



with the subscript 1 indicated that the Hamiltonian ap- 
plies to a single qubit, and the corresponding unitary 
operator is transformation 



Ui(t) — exp 
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We can now use Í7i and U2 to generate an entangled 
state. For example, the Bell state 
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is generated from the product state | J.)i ® | i) 2 by 



= e - ï7 r/4 eï (7r/4)<T x ( 2 ' e i(7r/4)<T z < 1 ) e -i(7r/4) C 
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In liquid-state NMR we do not begin with initial pure 
states, but we begin with a mixed state density operator 



(7) 



with H the system (individual nuclear spins on a sin- 

tmolecule) Hamiltonian, — (fegl 1-1 and Z = trp e 
. At high temperature, as is the case for present- 
day liquid-state NMR quantum computation |l5| , |24| , the 
state is very close to the identity so 
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where d = 2 N is the dimension of the Hilbert space for N 
qubits, 1 is the identity operator in this tensor product 
space and p\ is an arbitrary density operator. For exam- 
ple, in the case of a molecule with two spins and scalar 
coupling we have, 



with J <C u>i, tü2, and thus 
11 
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(9) 



(10) 



where p is the ratio of the Larmor freqüències of spin 2 
and 1. By a carefully tailored sequence of RF pulses, 
any two spin unitary transformation of this state can 
be achieved. Furthermore, using a spatially non-uniform 
magnètic field pulse (a gradient pulse) and averaging the 
varying phases over the sample we can effect particular 
non-unitary transformations j25j ]. With these two tech- 
niques it is possible to prepare the system in a so called 
pseudo-pure state of the form 
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1 + el*) (|*| 



(11) 



for |*)( V I'| a pure, and possibly entangled, state for 
N spins. It is possible to place bounds on the value 
of e for which the total state in Eq. (10) is entangled 
[ ]ï(| , that is a state which cannot be written as a convex 
combination of factorisable density operators. In typi- 
cal experiments e ~ 10~ 5 , a value which is too small 
for these states to be entangled. The pseudo-pure states 
produced in two qubit NMR quantum information pro- 
cessing experiments are not entangled and thus the spin- 
spin correlations at a fixed time have a purely classical 
interpretation. 
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Even though entangled states have not yet been pro- 
duced in NMR quantum information processing experi- 
ments, this does not mean that the system is not quan- 
tum mechanical. The important question is whcther 
there is a classical description of the dynamics of these 
experiments. This is a question that can and must be 
answered in a way that does not depend on the initial 
and final states of the system. It is a question concern- 
ing the propagator, or Greens function, for the dynamics, 
not the initial and final states. It is possible that the ini- 
tial and final states may exhibit no quantum correlations 
and have a perfectly vàlid classical description, yet not 
be connected by a classical dynamical model, be it deter- 
ministic or stochastic. In the case of stochastic dynamics 
the answer to this question will involve a specification of 
transition probabilities. As we show below, via a par- 
ticular but well motivated classical model, it is possible 
that the initial and final states are described by separable 
states (and thus could be interpreted as a perfectly vàlid 
classical probability distribution) yet no positive transi- 
tion probability exist to connect them either globally or 
infinitesimally. 

The quantum dynamics of a single spin-half system, 
described by the Hamiltonian Eq. (||), is equivalent to the 
classical dynamics of the corresponding magnètic dipole 
in an applied field |22|]l6|l . In both cases we find a linear 
precession of the average magnètic dipole about the di- 
rection of the applied field. Of course the states involved 
in the classical and quantum case are quite different. In 
the quantum case we can express the time-dependent 
state of the spin as 

p(t) = ~(l+n(t)-a) (12) 

with ft(t) a unit vector and a = e x à x + e y à y + e z à z . 
If the quantum Hamiltonian is H = B ■ 3 the quantum 
dynamics is given by the solution to the Bloch cquations: 

dn -> , 
~dt = ~ X (13) 

This is the same equation of motion as a classical (unit) 
magnètic dipole, n, in a magnètic field B. The equation 
describes the linear precession of a point, n, on the unit 
sphere around the direction of B at a rate \B\. 

Instead of a single classical dipole suppose we had a 
distribution of dipoles described by some initial proba- 
bility distribution function on the sphere, Qo(n). As the 
precession on the sphere is linear, each vector will precess 
at the same rate \B\ around the direction of B. The dis- 
tribution at time is then simply Qo(n(t)). In other words 
the solution of Eq. ( |Ï3|) are the characteristic equations 
for the equation of motion of the distribution function. 
The distribution simply rotates without distortion at a 
constant rate around B. 

However, the classical and quantum dynamics that re- 
sult for two magnètic dipoles interacting via the spin- 
spin interaction Eq. (Q), an entangling interaction, are 



very different as we show below. Thus we conclude that, 
while at present liquid-state NMR may not have access 
to entangled quantum states, it does allow us to realize a 
quantum dynamics for those states that is not be realized 
classically. It is the dynamics that is quantum in liquid - 
state NMR not the states. Liquid state NMR allows us 
to experimentally study the quantum dynamics of many 
coupled qubits and at present probably the most inter- 
esting element is to understand the amount of control 
we have on this dynamics. The corresponding classical 
system, with the same Hamiltonian, could never achieve 
this. 

To explain this we first note an equivalence between 
the spin-spin interaction of Eq. ([]]) and a nonlinear top 
model. Consider the collective angular momentum oper- 
ator S z defined by 

S z = \{^+*?\ (14) 

It is then easy to see that 

Sl = \(l + ^*^) . (15) 

Thus up to an irrelevant additive constant we obtain 

H 2 = J -Sl (16) 

The question of the dynamics is reduced to studying the 
quantum and classical dynamics of this nonlinear top. 
Note that this Hamiltonian commutes with S 2 the total 
angular momentum operator. Thus the system cannot 
evolve out of the subspaces corresponding to the irre- 
ducible representations of a two spin system. There are 
two such subspaces, the triplet with s — 1 and the sin- 
glet with s = 0. If we begin with the state in which both 
spins are down we cannot leave the triplet subspace with 
this Hamiltonian. Of course combinations of the two- 
spin and single-spin unitary operators will mix the two 
irreducible subspaces. However, the quantum and classi- 
cal dynamics that results from the nonlinear top interac- 
tion (with the same hamiltonian form), equivalently the 
two-spin interaction, are different regardless of the initial 
states as we now show. To be specific we will consider 
the dynamics restricted to the triplet (s = 1) subspace. 

We will follow closely the presentation of Sanders pfij ] 
concerning the classical and quantum dynamics of non- 
linear tops. We will assume that the physical interac- 
tion between two spins is fixed as the scalar coupling of 
two magnètic dipoles. The corresponding Hamiltonian 
is then fixed and we can compare the dynamics of ob- 
servable quantities that results when the interaction is 
treated either quantum mechanically or classically. It is 
conceivable that the exact quantum dynamics could be 
simulated exactly by a different classical Hamiltonian. 
After all we could always simulate the quantum system 
on a classical computer which is indeed a classical system 
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with a very complex time-dependent Hamiltonian. How- 
ever, we believe it is unlikely that any classical Hamilto- 
nian no matter how complex could simulate the quantum 
dynamics over a fixed time interval. We will return to tiris 
point in the discussion section below. 

The classical dynamics of a nonlinear top is deíined by 
the Hamiltonian 



H = LüS z 



J_ 

2s" 



S 2 



(17) 



where we have included a linear precession term with u> 
the linear precession frequeney. In this case the quan- 
tity S z is the z-component of the classical angular mo- 
mentum of the top. The first term describes the linear 
precession of the angular momentum vector about the z 
axis at the constant rate ui. The second term describes 
a nonlinear precession about the z axis at a frequeney 
that depends on the z component of angular momentum. 
The classical mechanics is described by the motion of a 
point in a spherical phase-space embedded in the thrcc- 
dimensional Euclidean space with coordinates S x , S y , S z 



with Si 



Sz — s 2 l 270. The classical states are 



probability distributions which describes an ensemble of 
tops with a distribution of angular momentum directions 
(every top in the ensemble has the same magnitude of 
total angular momentum). The points on the sphere of 
radius s are conveniently parameterized in polar coordi- 
nates as 



S/s — (sin 9 cos (f>, sin 9 sin </>, cos ( 



(18) 



However, we will use the stereo-graphic projection of the 
sphere onto the complex plane deíined by the map 



= e # tané)/2 



(19) 



The north pole (S z = s) is mapped to the origin (z = 0) 
and the south pole is mapped to infinity. The equator 
is mapped to the unit circle. In this conformal mapping, 
distributions with circular contours are mapped to distri- 
butions with circular contours in the complex plane. The 
dynamics of a distribution of points is easily described. A 
linear precession about the z axis simply causes a distri- 
bution of point to rotate about the origin in the complex 
plane without changing it shape. However the nonlinear 
precession causes the distribution to shear as different 
parts of the distribution with different vàlues for S z may 
have different precession rates. In the long time limit the 
distribution will tend to become smeared around the ori- 
gin in the complex plane (see Sanders |pr| for a pictorial 
representation). As we will see this is very different from 
what happens in the quantum case where the shearing 
ceases after some time and revivals and fractional revivals 
of the initial state oceur. 

In order to make a comparison with the quantum dy- 
namics we need to consider the dynamics of a distribution 
of points on the sphere. This is because a quantum state 
cannot be perfectly localized at a point on the sphere. 



The classical state of the system is described by a proba- 
bility distribution Q(z) of the vectors S corresponding to 
z. The expectation vàlues for the components of angular 
momentum are given by 



£{S Z ) 



d/j,(z)Q(z) 
dfi(z)Q(z 



*(i-M 2 ) 
i + M 2 

2sz* 



(20) 
(21) 



l + |z| 2 

where the integration measure in the stereo-graphic plane 



d M («) = ?Í±l(l + |z| 2 )- a 



(22) 



We have chosen the prefactor 2s + 1 as a sealing of the 
classical probability distribution, which makes the com- 
parison with the quantum case more convenient. 

The classical dynamics is described by a Liouville equa- 
tion 



dQ -{R r>\ 
-to-{H,Q} 



(23) 



where the Poisson bracket { , } may be determined us- 
ing {Si,Sj} — J2k e ijkSk- The Liouville equation is a 
first-order partial differential equation of the form 



dQ(z,t) 
dt 



d 

-z —Q(z,t) + c.c 
oz 



where the equations of motion are 
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i + M 
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(25) 



The solution to this is easily found after noting that \z\ 2 
is a constant of motion. Thus 



z{i) = exp 



it [ UI + J 



i-M 2 



i + M 2 

which in polar coordinates becomes, 
0(t) = 9 

(f){t) = (f) — Ult — Jt COS #o 



z(0) 



(26) 



(27) 
(28) 



where 0o ; $o are the initial vàlues. In this form it is 
particularly easy to see that the dynamics is a rotational 
shear of the sphere around the z axis. 

The solution for the probability density is 



Q(z,t)= d 2 z'T(z,z';t)Q(z',0) 



where the propagator is deíined by 

T(z,z';t) = S 2 (z'(t)-z) 



(29) 



(30) 



where z'(t) is the solution to the equation of motion for 
time t starting with the initial point z' . Whcn J = 
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we recover the previous result for a classical magnètic 
dipole: the distribution simply rotates, without distor- 
tion, at a constant rate ui about the z axis. The effect of 
the nonlinear term proportional to J causes a rotational 
shearing of the distribution around the z axis. 

We may include additional stochastic dynamics on top 
of the Hamiltonian dynamics. However, it is important 
to note that if Q{z, i) is a probability distribution then 
the propagators T must be positive and may be inter- 
preted as transition probabilities. The Hamiltonian Li- 
ouville evolution is a special case. As the propagator is 
simply the Greens function for the evolution equation the 
positivity requirement for the propagator restricts the al- 
lowcd form of evolution equations. It is well known that 
the allowed forms correspond to Fokker-Planck equations 
and can contain at most second-order derivatives with 
positive definite diffusion matrices |^8|] . In other words, if 
the propagators are to be positive the evolution equation 
is necessarily restricted regardless of the initial or final 
conditions. We may thus define allowed classical dynam- 
ics either in terms of positive transition probabilities or 
in terms of the differential operator for the dynamics. 

To compare the quantum and classical dynamics we 
now need to define a relevant quantum distribution. It is 
argued in references [ p6|]2S| ] that the appropriate object 
is the matrix elements of the quantum density operator 
in a coherent state basis. In the case of the harmònic os- 
cillator, these are the coherent states of the Heisenberg- 
Weyl group, and the resulting distribution is a true (ie 
positive) probability distribution for simultaneous mea- 
surement of position and momentum [ p0[ . In the case of 
angular momentum we can use the SU(2) coherent states 
pï[ defined by 



\z)=R{z)\s,s) l 



= (i + izi 2 r s £ 



where |s,m) are the 2s 
rotation operator is 



2s /2sN l/2 



z \s, s — m) 



(31) 
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1 eigenstates of S z , and the 



R{z) = exp — iOn.S 



(33) 



with the unit vector n = (sin0, cos <fi, 0). The states \z) 
are product states in terms of the qubits which are ro- 
tated from the state |0)i|0)2 by the angle cf> and 9 on the 
Block sphere: 



\z) 



1 



(i + M 2 ) 

The function 



(|0)i + *|l)i)®(|0) 2 + z|l) a ) 



Q(z,t)=tr(p(t)\z){z\) 



(34) 



is a true (that is positive) probability distribution for 
measurements defined by the projection operator valued 
measure (POVM) \z) (z\d 2 z/2n. Note that all allowed 



distributions are necessarily positive (and bounded) from 
the construction of Q(z, i) as a trace of the product of a 
positive operator and a projection operator. For exam- 
ple, the Q function for a particular atòmic coherent state 
\z ) is 



Q(z) 



(l + z^z)(í + z z*) 
(1 + z z*){l + zz*) 



(35) 



The first moments are given by integrals over the Q- 
function as 



(S x - iS y ) = / dfx(z)Q(z) 



2(s+l)z* 

' TTW 



(36) 



d^z)Q(z) ( {S+ l ) + ] J Z?ï ) (37) 



The important point to notice is that the averages of 
S± are given by the same functional form as the clas- 
sical case, apart from an additional term that becomes 
negligible in the semi-classical limit as s — > oo. The 
second-order moments are given by 



<£ 2 > 
(S 2 Z ) 



d[i(z)Q(z) 



2(2s + 3)(s + l)z* 



(i + M 2 ) 2 

d/j,(z)Q(z) 

{s + l) 2 - 2s(s + 2)\z\ 2 + (s 2 + 4j + 5)|z| 4 



(l + \z\ 2 ) 2 



Thus, even though the Q function is a true probability 
distribution, its marginals do not give the quantum ex- 
pectation vàlues: an additional rule is needed to connect 
averages over the Q function to the quantum averages. 
This is analogous to the case for the harmònic oscilla- 
tor coherent states [^9| . I n the case of a spherical phase 
space, however, the difference appears already at the level 
of the first-order moments. 

Taking matrix elements of the quantum Liouville equa- 
tion 



we obtain the evolution equation 



(40) 



d_ 

dt 



Q{z,t) 



M 2 



i + * 



j_ d_ 

2~s Z ~d2 



xz—Qíz, t) + c.c 
oz 



(41) 



This equation is linear in Q; thus, p\ in Eq. (jg) will obey 
exactly the same equation. In the limit of s — > oo, with 
ui, A held constant, the equation reduces to the first-order 
differential equation of classical dynamics. The differ- 
ence between quantum and classical dynamics is due to 
the second order differential operators. Note that while 
these terms are second order, they are certainly not of 
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the kind expected for a diffusion equation, as the corre- 
sponding diffusion matrix would not be positive defmitc 
This is a familiar feature of the difference between classi- 
cal and quantum dynamics as reflected in the dynamics 
of a quasi probability distribution and was first noted in 
the context of quantum òptics |52| 

In rcfercncc p6| it is was shown how extreme this dif- 
ference could be. For example, at times t = irs/J an 
initial coherent state (or coherent pseudo-pure state) \ Z q) 
would evolve, in the rotating frame (u> = 0) into the pure 
state 



(46) 



2- 1/2 (e- M / 4 |zo) + (-l) s e"/ 4 |-z )) 



(42) 



for which the resulting Q-function is double peaked. The 
state is entangled in terms of the pure state representa- 
tion. No state such as this could ever be obtained from 
the classical dynamics given in Eq. ( |29| ) starting from the 
initial state Eq. (|35|). 



III. DISCUSSION 

We noted above that the difference between the quan- 
tum and classical dynamics appears through the second- 
order derivatives in Eq. ([ïï|). How does this difference 
become manifest in the propagator for the equation? It 
is far from clear that we can find a positive propagator 
corresponding to the differential operator on the right 
hand side of Eq. pi]), especially as this equation is not 
of Fokkcr-Planck form (the diffusion matrix is not posi- 
tive definite) . Nonetheless all initial and final Q-functions 
must be positive by construction. We now show that the 
quantum evolution of the Q-function cannot be written 
in terms of transition probability propagators. Despite 
this it can be obtained uniquely from the initial condi- 
tion, Q(z, t = 0) 

To begin we define the Q-function amplitude for a pure 
state \tp(t)) 



Q(z,t) = (z\i>(t)) 



(43) 



The Q-function is given by the modulus squared of the 
Q-function amplitude, Q(z,t) = \ Q(z,t)\ 2 . The linearity 
of the Schròdinger equation now requires that the dy- 
namics of the Q-amplitude can be written in terms of a 
linear propagator 



Q(z,t) = / d 2 z'C{z,z';t)Q{z',0) 



(44) 



where Q(z,0) is the Q-function amplitude for the initial 
state and the propagator (for s = 1) is constructed from 
the unitary time evolution operator U(t) as 

C(z, Zl ;t) = (z\U{t)\ Zl ) 
1 



(l + |*| a )(l + N 2 ) 




Thus the Q-function dynamics can be written in integral 
form as 

Q(z,t) = J d 2 z 1 d 2 z 2 £(z,z 1 ;t)C*(z,z 2 ;t) 

xQ(zi,0)*Q(z 2 ,0) (47) 

In the case of a general initial state p(0) the equation 
becomes 

Q(z,t) = í d 2 z 1 d 2 z 2 C{z,z 1 -t)C*(z 1 z 2 -t){ziW)\z2) 



(48) 

This last expression seems to suggest we need to know 
more than just the initial condition Q(z, 0), but this is 
not the case. The matrix elements of p in the coherent 
state basis suffice to uniquely determine the state and 
thus uniquely determine the off-diagonal matrix elements 
pï| , by analytic continuation. A similar statement may 
be made about the propagator in Eq. (^|) . We only need 
to know the diagonal matrix element of U^\z)(z\U in or- 
der to determine the total propagator in Eq. (fí8|). Thus 
knowledge of the initial Q-function Q{z, 0) and a positive 
linear propagator K.(z,Zi;t) — (zi\U(t)^\z}(z\U(t)\z 1 ) 
uniquely determine the solution to the quantum evolu- 
tion equation for the Q-function. This is of course also 
true for the classical evolution equation, through a sim- 
ple convolution of the initial state and the propagator. 
However, the propagation integral in the quantum case, 
Eq. (^8|), is not a positive (or even real) function and can 
have no interpretation as a transition probability. 

It is generally accepted that uncontrollcd interactions 
with an environment enable the quantum and classical 
dynamics to be reconciled when states of the environ- 
ment are averaged over [^3) . In the case of two coupled 
spins, there are a variety of possible environmental inter- 
actions. In order to illustrate the principle of decoher- 
ence, we take the simplest possible case in which both 
spins are coupled equally to the same environment with 
a Hamiltonian that conserves the total z-component of 
angular momentum, S z . While this collective dephasing 
model is not very realistic for NMR experiments it will 
illustrate how decoherence can cause the quantum prop- 
agator in Eq. ([Ï8]) to become diagonal. 

The collective dephasing master equation is given by 

ü 

p = Vp = -i[H,p]-^[S z ,[S z ,p}} (49) 
Zs 

The general solution for the Q-function may then be writ- 
ten as 



Q(z,t) 



d 2 Zl d 2 z 2 V(z;z 1 ,z 2 ,t)(z 1 \p(0)\z 2 ) (50) 



7 



where the propagator is given in terms of the coherent 
state matrix elements of the dynamically propagated off- 
diagonal projector, 



V(z;z 1 ,z 2 ,t) = {z\e m (\z 1 ){z 2 \)\z) 



(51) 



For short times we can expand this to linear order in t. 
The dominant non-Hamiltonian terms in total spin s are 
given by 



V{Z\ Zl,Z2,t) 



1 



7SÍ 



Z2 



2 V(l + |zi| 2 )(l + |z 2 | 2 ) 
x'P(z;zi,z 2 ,0) + ··· 



(52) 



The cofactor of jst/2 is less than unity. For two spin half 
systems with s = 1 the coherence decay rate is then set 
cntircly by 7 which could be small compared to the time 
scale set by the cohernet interaction, J _1 . When the ini- 
tial state is concentrated on \zi\ = \z 2 \ the decoherence is 
also small. The coherence decay rate for large semiclas- 
sical systems for which s » 1 can be large. Of course 
the resulting propagator is not the same as the classical 
Hamiltonian result (Eq. (|29|)) in the absence of dephas- 
ing as the coupling to the environment would add some 
level of phase diffusion to the classical dynamics as well, 
broadening the delta function in Eq. (p9|). In the case 
of s = 1 considered here one would need to consider all 
the terms in the short time expansion for V(z; zi, z 2 , í), 
but it remains the case that terms with large vàlues of 
\zi \ — \z 2 \ are rapidly suppressed. 

Other two-qubit interactions might be considered, for 
example the exchange interaction |p4|, 



H^= J -a^.a^= J -0.S-\) 



(53) 



where the total spin operator is S = (a^ + a^)/2. This 
is the Hamiltonian for a rigid body wiht an isotropic 
moment of intertia with no external torques and with 
moment of inertia J . The classical phase space is a 
cylinder with canonical coordinates S on the axis of the 
cylinder and 9 around the circumfcrcncc. The classical 
equations are the same as those of a free mass with pe- 
riòdic boundary conditions in position at ±7r : 9 = Js 
and 5 = 0. The first of these again indicates a rota- 
tional shearing of a localized distribution on the phase 
space. When the shearing causes different parts of the 
state to overlap on the classical phase space, we expect 
the corresponding quantum system to exhibit interfer- 
ence fringes. Thus the quantum and classical dynamics 
of two interacting magnètic dipoles with such an interac- 
tion must differ. (Clearly such an interaction would need 
a different mechanism in classical systems as exchange 
interactions are quantum mechanical.) 

As we have seen the nonlinear top can generate the su- 
perposition state (e -I7r / 4 |s, s) z + e~ l7r / 4 |s, — s) z ). A simi- 
lar state can be generated by a sequence of CNOT gates 
on N — log 2 (2j + 1) qubits. A product state of N qubits 



can be written in terms of a binary string \X) = Iíf\xi) 1 
where Xi is the i'th term of the string X. Alternatively 
X could encode an integer k in binary form. The maxi- 
mally entangled state of N qubits, |000 ... 0) + 1 111 ... 1) 
would then take the form |0) + \M) where M = 2 N . 
Such as state is easily generated in a quantum computer 
by a single Hadamard gate on the first qubit followed 
by a cascaded sequence of controlled NOT gates. If we 
change the notation for angular momentum states so that 
\s, m — s) z = \m) where m = 0,1,... 2s, then the an- 
gular momentum superposition state generated by the 
nonlinear top is equivalent to the maximally entangled 
state. This equivalence suggests that a nonlinear top 
may itself be made to act as a quantum computer in 
the 2s + 1 dimensional Hilbert space. Indeed numerical 
evidence exists p5|] that a time-dependent Hamiltonian 
of the form H(t) = B(t) ■ S + J(t)S^ can generate any 
state in the Hilbert space by a suitable choice of the time 
dependent coefhcients. We conjecture that a sequence 
of pulses U(6 n ,<p n ,Xn) = exp(-i6>„n„ • S) exp(-zx„5' 2 ) 
where n = (cos <f>, sin</>, 0) will suffice to perform quan- 
tum algorithms in the Hilbert space of the nonlinear top. 
Further work remains to be done to determine whether 
there are interesting efncient algorithms in terms of num- 
ber of pulses and (say) log(s). 



IV. CONCLUSION 

We have investigated classical and quantum Hamilto- 
nians of the non linear top and show that they producc 
different observables even in the presence of highly mixed 
states. This shows that even if a state is separable its 
classical and quantum evolution can differ. Thus this for- 
tifies the claim that the evolution is one of the quantum 
elements of liquid-state NMR QIP. However it does not 
prové that bulk ensmeble NMR quantum computation is 
uniquely quantum. We have used a particular classical 
model, that is one motivated by the physics of two in- 
teracting magnètic dipoles, to compare the quantum and 
classical dynamics. It may well be the case that for states 
that are close to the maximally mixed state there is an- 
other classical model (in effect a hidden variable model) 
that correctly describes the dynamics as far as is requircd 
to model the observed results in NMR experiments ]Ï6|| . 

Even if the transformation of the state is not classical 
our present work does not show that using these mixed 
state qubits is as powerful as a pure state quantum com- 
puter. If we use pseudo-pure state as has been done in 
present experiments, the signal to noise decreases expo- 
nentially, thus rendering the algorithms inefficient. In 
this case liquid-state NMR does not provide any advan- 
tage (with respect to speed) over classical computers. In 
the absence of noise, Schulman and Vazirani |3l| have 
shown how to efncicntly transform the initial mixed state 
into a pure state. Some simulations of quantum systems 
evolving under unitary transform and the algorithm sug- 
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gested in |}7| could also be implemented efficiently as 
long as noise is negligible. Thus liquid-state NMR (in 
the absence of noise, i.e. on time scales smaller than 
T2) offers a testing ground for QIP. The dynamics is the 
correct one. Moreover it has the right error model to pro- 
vide a test bed for small quantum computations (up to 
roughly 10 qubits). The difference between liquid state 
NMR and other proposals for quantum computation is 
that, in the former case we do not have, even in prin- 
ciple, a method to make the quantum evolution robust. 
In practice this does not distinguish NMR from other 
present devices as none have approachcd the accuracy 
needed for scalability. 

In the present work we have neglcctcd the possibil- 
ity of hidden variables. The reason is that there is a 
trivial hidden variable model which would explain all the 
present experiments in QIP (including not only NMR but 
also other technologies) as long as observations are not 
made on space-like surfaces. The model can be thought 
of as a classical computer which simulates the quantum 
evolution and telis the bits of the physical system how 
to behave so they mimic quantum mechanics. Although 
this model can describe all the experiments in QIP today, 
the amount of resources it uses compared to its quantum 
counterpart always seem to be exponential in the size 
of the problem (the model of Schack and Caves is such a 
model) . It is important that the resources account for not 
only the signal to noise but all resources as it is usually 
easy to trade one resource for another. We would like to 
conjecture that what distinguishes classical and quantum 



computation is the amount of algorithmic information re- 
quired to produce the answer of some problem starting 
from a fiducial state. Classical devices require an expo- 
nential amount of algorithmic information compared to a 
quantum computer to answer certain problems. Simulat- 
ing liquid-state NMR experiments is one such problem 
compared to known classical algorithms. 

So where did the power of quantum computation comc 
from? In this paper we have argued that the power comes 
from the dynamics of the system. A similar argument 
has been made by others JÏ(J,[Ï6|]. We have criticized the 
view that entanglement is the source of power of quantum 
computation by giving algorithms and dynamical evolu- 
tion which do not depend on entanglement. There may 
be a variety of elements which make these devices more 
powerful, and for another unusual method to quantum 
compute see However, as long as we lack a proof 

that we cannot simulate quantum systems efficiently, it 
is hard to attribute the source of the power of quantum- 
ness, and we await more powerful arguments. 
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(z+z*)/2 




(z+z*)/2 

FIG. 1. Comparison between the classical and quantum evolution of the Q(z) function. The initial state is given by a). It 
corresponds to the distribution for the quantum state \z = 1). b) and c) depict the classical and quantum distribution at time 
t = 2tï/.J (for uj = 0). The classical evolution follows the equation of motion Eq. ( pi| ) and the quantum one Eq. @. The 
discrepancy between the classical evolution and the quantum one is evident. 
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